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ABSTRACT 


Parameterizations of the frontal area index and canopy area index of natural or 
randomly distributed plants are developed, and applied to the estimation of local 
aerodynamic roughness using satellite imagery. The formulas are expressed in terms of the 
subpixel fractional vegetation cover and one non-dimensional geometric parameter that 
characterizes the plant's shape. Geometrically similar plants and Poisson distributed plant 
centers are assumed. An appropriate averaging technique to extend satellite pixel-scale 
estimates to larger scales is provided. 

The parameterization is applied to the estimation of aerodynamic roughness using 
satellite imagery for a 2.3 km 2 coniferous portion of the Landes Forest near Lubbon, 
France, during the 1986 HAPEX-Mobilhy Experiment. The canopy area index is estimated 
first for each pixel in the scene based on previous estimates of fractional cover obtained 
using Landsat Thematic Mapper imagery. Next, the results are incorporated into Raupach's 
(1992, 1994) analytical formulas for momentum roughness and zero-plane displacement 
height. The estimates compare reasonably well to reference values determined from 
measurements taken during the experiment and to published literature values. The approach 
offers the potential for estimating regionally variable, vegetation aerodynamic roughness 
lengths over natural regions using satellite imagery when there exists only limited 
knowledge of the vegetated surface. 


Keywords: Canopy area index; Aerodynamic roughness; Remote sensing; Coniferous forest 



1. Introduction 


Most model parameterizations of the surface momentum, sensible heat and latent heat 
fluxes, and many field and remote sensing techniques to measure them, rely on Monin- 
Obukhov similarity theory (e.g., Brutsaert, 1982; Stull, 1988). In the case of momentum, 
the similarity equation is expressed, 


u(z) 


u* k‘ 


In 


z - d„ 


- v(C) 


( 1 ) 


where u(z) is the surface layer wind speed at height z, u* is the friction velocity, k is von 
Karman's constant, and t| r(t) is the adjustment function due to stability. The momentum 
roughness length, z 0 , is conceived as the zero velocity intercept of the vertical mean wind 
profile. The zero plane displacement height, do, is an elevation adjustment due to the 
conceptual vertical shift of the logarithmic profile in the presence of the roughness elements. 
Analogous similarity expressions are used for the sensible heat flux in terms of z 0 h, the 
roughness length for sensible heat, and the latent heat flux in terms of zov, the roughness 
length for water vapor. 


Numerical atmospheric simulation models used today for short-term operational weather 
forecasting including mesoscale models, or long-term climate analysis including general 
circulation models, employ some form of the above similarity equations. In order to run the 
models, roughness lengths are prescribed for each grid box in the domain prior to the 
simulation. However, since roughness data are not readily available at that resolution, 
values are assigned by associating roughness with vegetation type or land classification 
using look-up tables (e.g., Grell et. al., 1994). In the case of most operational general 
circulation models, Zq, Zoh, and Zqv are assigned the same value despite evidence to the 
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contrary (e.g. Brutsaert and Sugita, 1992). The method usually involves the use of climatic 
land classification values that vary only on a seasonal basis but do not change from year to 
year. 

The reason for the above approach lies in the difficulty in estimating roughness. 
Theoretically, z Q and do are independent of flow and only a function of the geometry of the 
surface (Brutsaert, 1982; Monteith and Unsworth, 1990). While analytical formulas and 
numerical studies have been presented based on drag expressions (e.g. Lettau, 1969; 

Mason, 1988; Seginer, 1974; Shaw and Pereira, 1982; Raupach, 1992, 1994, Wooding et. 
al., 1973; see summary by Brutsaert, 1982), they require knowledge of mean surface 
characteristics such as height, frontal area, surface density, and surface geometry. For 
vegetated surfaces in particular, key variables are the mean canopy height, frontal surface 
area and the mean plant density. Unfortunately, those characteristics are generally 
unknown, especially for natural regions where plants are randomly dispersed, that precludes 
application of the above formulas. 

Practically, zq and do for various land surface geometries are best determined from site- 
specific field experiments from the analysis of flux-profile measurements during neutral 
conditions, a method that requires substantial investment in manpower, time, and 
instrumentation. Over time, values have been tabulated for a wide range of representative 
surfaces (e.g., Brutsaert, 1982; Stull, 1988). 

Improved characterization of the land surface, using satellite imagery, may lead to 
improved estimates of aerodynamic roughness lengths and significant progress over present 
applications. Over the past two decades, satellite remotely sensed observations of reflected 
visible and near-infrared solar energy have been used to estimate an increasingly greater 
number of plant attributes such as leaf area index (LAI), biomass, and fractional cover. 
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Most techniques involve correlation with empirical indices constructed from multispectral 
bands (e.g. Perry and Lautenschlager, 1984; Tucker ef. al., 1983; Asrar et. al., 1984; 

Huete, 1988), although sensitivity studies have shown their limitations (Jasinski, 1990; 
Huete and Liu, 1994). Other investigators have modeled the radiation physics through the 
plant canopy (e.g. Ross, 1981 ; Myneni et al ., 1989; also see summary in Asrar, 1989). 
However, that approach is often constrained by the relatively large number of model 
parameters, especially in heterogeneous regions. 

The inversion of geometric canopy reflectance models using satellite multispectral 
imagery has provided a practical approach for investigating regional vegetation variability, 
and a possible avenue for characterizing roughness. Plants are conceived as randomly 
distributed, three-dimensional elements superposed on a flat surface (Li and Strahler, 1985), 
with only a few key bulk plant parameters, often in non-dimensional form (Jasinski and 
Eagleson, 1990). The approach has successfully provided estimates of fractional cover and 
LAI over forested and semiforested landscapes using satellite imagery (Li and Strahler, 

1985; Franklin and Strahler, 1988; Jasinski, 1996). It is especially applicable in natural 
regions where the canopy cover often is spatially random. Since the method relies on the 
characterization of the bulk canopy geometry, it offers an opportunity for estimating the 
local aerodynamic roughness that also is a function of those same physical properties. 

2. Background 

2.1 Parameterization of zo and d 0 in terms of plant frontal area density 

One analytical approach for parameterizing the local momentum roughness length in 
terms of geometric canopy parameters is that by Raupach (1992, 1994). Using dimensional 
analysis and physical hypotheses on the partitioning of surface stress, he derived an 
expression for z<yh in terms of the drag on the roughness elements and on the underlying 
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surface. The parameterization provides "a relatively simple unified theory" (Mahrt, 1996) 
for turbulent transport within canopies. 

A key plant parameter is the frontal area index or density, X, defined as the ratio of 
frontal area of roughness elements, from the mean wind direction, per unit ground area. 
For solid elements. 


_ nbh 

S (2, 

where bh represents the frontal area (breadth x height) of a single roughness element, n is 
the number of elements, and S is the total ground area. Raupach (1994) later determined it 
convenient to replace X with the canopy area density or index. A, or the total single sided 
area of all the canopy elements per unit ground area. It differs from the LAI in that the latter 
usually includes only the green leaf area, while A includes all canopy elements that absorb 
momentum. For isotropically oriented elements, the relationship between the two is 
approximately A = 2X (Raupach, 1 994). 


Raupach's (1994) aerodynamic roughness equations are expressed 


( <V 


f u h „ 1 
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where Uh is the velocity at the top of the canopy, Wh is a velocity profile adjustment, and 
the remaining terms are as defined above. The ratio of u*/Uh is expressed 
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where C s and Cr are the drag coefficients for the ground and roughness elements, 
respectively, c is an empirical coefficient related to wake spreading, and (u*/Uh)max is the 
maximum observed ratio when the flow begins to skim over the top of the canopy instead of 
penetrating to the surface. Raupach (1994) developed an empirical expression for the 
displacement height , or 


\_ 

d 0 1 -exp [-(c dl A)^J 
h l 

(c d iA ) 2 

where Cdi is derived from a wide range of laboratory and field data. The equation is 
designed to approach the proper limits as A approaches zero and infinity. 


2.2 Geometric canopy model for natural or randomly distributed plants 

The geometric model, shown in Figure 1 , is characterized by both its shape and its 
spatial distribution. The canopy itself is represented as three-dimensional, geometric 
elements positioned on a flat horizontal surface. The plants are assumed to be geometrically 
similar, that is, the individual canopies within a plant stand, although of different size, are 
considered to be one geometric shape with a constant width to height ratio, D/H. Various 
common shapes such as cones, cylinders and ellipsoids can be employed as shown in 
Table 1. 
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Spatially, the model assumes a random horizontal plant distribution. Plant centers of 
natural regions are most often described by the Poisson process (Diggle, 1983; Whittaker, 
1975). Further, plants with large diameters or situated in dense stands tend to overlap one 
another. Consequently, the mean vertically projected canopy area per unit ground area, or 
the fractional canopy cover, m, is mathematically expressed after the probability work by 
Kellerer (1983), or 


m = 1 - exp(-Atp) (6) 

where At is the mean vertically projected canopy area of a single tree, and p is the spatial 
Poisson density or. 


p = n/Ap. (7) 

where Ap is the unit area, in this case the area of the satellite pixel. 

A final characteristic of the geometric model, necessary for remote sensing, is the 
parameterization of subpixel fractional ground shadow cast by the canopies in terms of their 
geometric shape, spatial density, and the solar zenith angle. It was previously shown 
(Jasinski, 1990) using Equation (6) that for any Poisson distributed plant, a 
parameterization of the subpixel ground shadow exists for satellite pixels greater than 
several tens of meters, or 


gs - l-m-(l -m)1 + l. 


( 8 ) 
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where m is the total fractional canopy cover in a given pixel, g s is the total fractional ground 
shadow in a pixel, and T| is a non-dimensional, solar-geometric similarity parameter defined 
as the ratio of the mean shadow area cast by a single tree, A$, to its mean projected canopy 
area or. 


*1 = 



(9) 


Using T| in this non-dimensional form allows one to generalize the shadow parameterization 
for any canopy geometry. Analytical expressions for T| for several geometrical shapes are 
shown in Table 1. The expressions are developed in terms of the solar zenith angle, 0, and 
the ratio of some maximum plant width to plant height, or D/H. Qualitatively, large rj 
represents tall, bulky plants that cast large shadows while small T) represents thin short 
plants. The T| can be estimated from Table 1 or directly from aerial photographs during 
clear-sky conditions. 


3. Parameterization of frontal and canopy area indices for natural regions 

3.1 Frontal area index for spatially random solid geometric canopies 

In order to apply Raupach's equations to satellite imagery, it is first necessary to 
develop a subpixel parameterization of the frontal area or canopy area index. Since the scale 
of the plant canopy (order of several meters) is much smaller than the scale of the typical 
satellite pixel (order of tens of meters or greater), variations in frontal density occur at the 
subpixel level. Consequently, a parameterization is required that relates frontal density in 
terms of a subpixel quantity that can be measured from satellites, such as fractional 
vegetation cover. Since the method will be applied mostly to natural or undisturbed regions. 
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it is also desirable to incorporate the spatially random nature of the plant canopy in the 
formulation. 

The present approach will follow the method used to parameterize subpixel shadows 
described above. The starting point is the equivalent definition of the frontal area density, 
X, or 


X = 



( 10 ) 


where Af is the mean frontal area of an individual plant or roughness element. The 
nondimensional frontal area similarity parameter, v, is now introduced as the ratio of the 
mean frontal area of an individual plant to the mean vertically projected canopy area or 

Af 


Physically, large v represents tall plants with narrow canopies while small v represents short 
plants with wide canopies. Theoretical expressions for v have been derived for various 
geometric shapes in terms of the geometric similarity ratio described earlier, as shown in 
Table 1. Combining Equations (10) and (11) yields, 


X = 



( 12 ) 


Assuming a Poisson horizontal distribution within each pixel, the combination of (6), (7) 
and (12) yields a subpixel parameterization for the frontal area index in terms of fractional 
cover and the non-dimensional quantity v, or 
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X = -v ln(l-m). 


(13) 


3.2 Canopy area index for spatially random solid geometric canopies 


The canopy area density is the total, single sided area of the canopy per unit ground 
area. For geometric solids, unlike leaf and branch elements, the assumption of isotropic 
components is not generally applicable. Therefore, the assumption that A - 2X can not 
exactly be applied except perhaps for certain (i.e. ellipsoidal) canopies. However, an 
expression for A can be formulated in a similar manner to the frontal area index. The 
canopy area similarity parameter is first defined. 



(14) 


where Ac is the total exposed (i.e. outside) surface area of the geometrically shaped canopy. 
For consistency with Raupach's (1994) formulation for momentum roughness, the "single- 
sided" canopy area is defined as one-half the exposed surface area or A c /2. Using the same 
approach as in Equations (10) through (13) in Section 3.1 above, the subpixel canopy area 
index parameterization becomes, 


N 

A = ln(l -m) 
2 


( 15 ) 
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Expressions for N for various geometric shapes are provided in Table 1. For canopies 
conceptualized as solid geometric elements, the relation between the canopy area and frontal 
area indices is, 


A = 



(16) 


where A = 2X. only when N/v = 4. 

Equations (13) or (15) provide the essential link between the estimation of fractional 
cover using remote sensing, and Raupach's method for estimating aerodynamic roughness. 
They indicate that if the canopy frontal similarity parameter can be prescribed in geometric 
terms of v or N, then once fractional cover is known from the satellite analysis, the canopy 
frontal density required for Raupach's equation can be estimated for each satellite pixel. 

A log-linear plot of the theoretical canopy area index versus fractional cover is provided 
in Figure 2. The figure is not specific to a particular canopy shape, but presented in non- 
dimensional form, is applicable to any plant with Poisson distributed centers. It can be 
utilized in conjunction with any plant shape when configured in the non-dimensional form 
using the canopy area similarity parameter, N. A similar graph could be plotted for X versus 
m. 


The physical interpretation of Figure 2 can be explained as follows: For sparse 
canopies, that is fractional covers less than about 0.20, the canopy area index increases 
roughly proportionally with increasing fractional cover, due to the behavior of ln( 1 -m) for 
small m. Physically, this is due to the fact that there is little overlap among the canopies at 
such low horizontal densities. However, for large fractional covers, it takes an increasingly 
greater number of trees for each additional unit of canopy area due to increasing overlap 
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among the dense trees. Consequently, the canopy area index increases less rapidly per unit 
increases in fractional cover at this higher range. 

The use of the Table 1 in conjunction with Figure 2 provides a mechanism, at least 
theoretically, to compare the canopy area indices of differently shaped roughness elements. 
For instance, roughness elements with a similarity shape of N = 2 and cover density 
m = 0.6 would theoretically possess the same canopy area index as roughness elements 
with a similarity shape of N = 10 and cover density m = 0.16. 

In some cases, previous knowledge from the shadow analysis can be used in the 
estimation of v or N. For example, knowledge of T) and the solar zenith angle can provide 
the ratio D/H in order to estimate v or N. 


4. Characterization of aerodynamic roughness of the Landes Forest 

4.1 Approach 

The estimation of aerodynamic roughness as presented herein is formulated as a three- 
step approach. First, the subpixel fractional canopy cover is estimated for each pixel in the 
scene. Second, the corresponding canopy area index for each pixel is estimated using the 
above formulas and fractional cover estimates. Third, the canopy area index is incorporated 
into Raupach's formulations in order to estimate roughness. Although other roughness 
formulas can be used, Raupach's seem the most straightforward. A simple averaging 
technique for the current problem is presented. 

As indicated in the Introduction, techniques for the first step in estimating fractional 
cover using geometric modeling are well documented and therefore will not be repeated 


13 



here. Thus, the present analysis will focus on the latter two steps where, building on 
previous work, they are applied to a region where published satellite-derived fractional 
cover estimates already exist. 

4.2 Estimation of subpixel fractional cover and leaf area 

One region that has received considerable attention is the Landes Forest near Lubbon, 
France, site of the June 1986 HAPEX-Mobilhy First ISLSCP Experiment (Andre et al, 
1986). This portion of the Landes Forest is flat and dominated principally by maritime pine 
(Pinus Pinaster Aiton) stands with an understory consisting of bracken fern. Interspersed 
between the pine stands are shrub and agricultural lands comprising less than one-quarter of 
the total area. 

Fractional cover estimates previously were made on a pixel by pixel basis at 30 m 
resolution for an approximately 2.3 km 2 portion of the region using nadir viewing Landsat 
Thematic Mapper imagery (Jasinski, 1996). The frequency distribution of fractional cover 
is shown in Figure 3. The mean fractional cover of the total area was estimated to be 0.67 
and the standard error was estimated at 0.056. 

The leaf area index of only the tree canopy, LAIc, without including the bracken 
understory, was estimated to be approximately 2.3 based on measurements at one tower site 
(Gash et al., 1989). Satellite-based measurements (Jasinski, 1996) for pixels within the 2.3 
km 2 Lubbon region possessing nearly complete canopy cover (m ~ 1) yielded a total leaf 
area index of both tree canopy and the bracken understory, or LAI, of about 3.1. Assuming 
the bracken typically possesses an LAI of about 1 .0, the satellite estimates of total LAI are 
consistent with the LAI C value observed by Gash et al. (1989). 
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4.3 Characterization of canopy geometric properties 


The canopy geometric properties required for the analysis were based only on the 
maritime pine, since it dominates over three-fourths of the region. The majority of the 
maritime pine stands in the Lubbon area were fairly mature, ranging in size with an average 
total height of about 20 m and an average density of approximately 430 stems ha 1 (Gash et 
al., 1989). Although no field measurements of the canopy geometry were made, analysis of 
aerial photographs indicate that the structure of the typical tree can be characterized as a 
bulky, conifer canopy situated on a long stem with a mean width, D, of approximately 5.1 
m. From a geometric modeling perspective, the trees are best represented as cones on a 
post, shown in Table 1 . 


In order to proceed, the ratio D/H, where H is the vertical height of the canopy itself, 
and the height of the stem below the canopy, h s , need to be estimated. Those values can not 
be determined directly from a single, nadir viewing aerial photograph. However, for this 
case their values can be estimated by knowledge of the shape of the tree and its shadow 
regime. Recalling previous estimates of t| = 1.05, and the solar zenith angle at time of 
overpass, 0 = 31.6° (Jasinski, 1996), and assuming geometrically similar trees, the 
definition of r| from Table 1 becomes. 


where 


and 


t) = 1 + 


n 


K 


cot X - 2 + x 


--{cos ‘(^-[nO-n^]} 
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.2x2 


X Sin ^Htane ) 
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(17) 
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Combining Equation (17) with the above limited information (i.e. D = 5.1 m, and H + h s = 
20 m) yields estimates of H = 12.8 m, D/H = 0.40, and h s = 7.2 m. A summary of all the 
Landes Forest site parameters is provided in Table 2. 

4.4 Estimation of canopy area index 

Continuing with the assumption of cones on a post, the value of N can be estimated 
using the formula in Table 1 , or, 

N = 1 + 

= 6.1 

Given N, Equation (15) is used to calculate the canopy area index for each pixel in the 
scene. The resulting frequency distribution over the entire scene is provided in Figure 4. 

The mean canopy area index is computed by first averaging over the fractional cover 
density in order to account for scaling effects, or 



where t is the number of pixels, and mj is the fractional cover of an individual pixel. For the 
present case, the results yield A= 3.4 for the entire 2.3 km 2 area. Following the same 
reasoning, it can be shown that the mean frontal area index is X.= 1.8 . Those values 
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seem reasonable and are consistent with analogous estimates for a pine forest (e.g. Garratt, 
1977). 

4.5 Application of Raupach's formulation to estimate of Zoth and doth for the Landes 
Forest 


The roughness quantities, zjh. and d<yh, theoretically can be estimated using A, 
Equations (3) through (5), and some knowledge of several constant parameters including 
Cr, Cs, (u*/Uh)max> c w , and cji. Raupach (1994) assumed global values for the constants 
based on field and wind tunnel data including forest sites, yielding the graphs in Figures 5 
and 6. Applying those same global values without any special fitting for each pixel in the 
Lubbon Landes Forest, together with the estimated canopy area index, yields the frequency 
distribution of Zc/h shown in Figure 7, and the frequency distribution of d<yh shown in 
Figure 8. 


In the present case, the forest area of 2.3 km 2 is relatively small, and h has been 
assumed constant. The mean roughness, z o^ can be approximated simply by 



h 



( 20 ) 


where fj is the frequency of z 0 j/h from Figure 7. Particularly in cases where h is not 
constant, a more sophisticated averaging procedure may be necessary. 

Applying Equation (20) to the frequency distribution in Figure 7 yields z c/h = 0.061 . 
This can be compared to the estimate obtained simply by applying the mean canopy area 
index (i.e. for A= 3.4 ) to Equations (3) and (5) that yields z c/h = 0.062. Thus, the 


17 



averaging approach is not significant for the present case. Applying either approach to the 
displacement height yields a mean ^o^ = 0.79. Thus, the results indicate a scene average 
of z o= 1.2 m and = 16.0 m. 

5. Discussion 

5.1 Comparison with reference values 

Roughness lengths and displacement heights have been estimated for the Landes forest 
using at least two independent field methods. Parlange and Brutsaert (1989) analyzed high- 
resolution radiosoundings under neutral stability over the region during the HAPEX- 
Mobilhy Experiment. They assumed that do is approximately five times as large as Zq [e.g., 
Brutsaert (1982), pages 113-1 16]. The analysis yielded z 0 =1.2 m and d o =6.0 m. On the 
basis of Paeschke's (1937) estimate that z 0 =h/7.5, they concluded that the effective height 
of the roughness elements felt by the atmosphere was near 10 m. This is intermediate 
between the average forest height of 20 m and the canopy heights in the clearings. During 
the same period. Gash et al. (1989) analyzed energy flux measurements made with eddy 
correlation equipment stationed on a tower 29 m above the ground. They estimated Zq to lie 
at a height between 1 .4 m and 2.4 m, or approximately 1 .9 m, and they estimated that do is 
about 75% of the mean tree height, or about 15 m. 

For canopy height taken to be the mean tree height of 20 m, the Parlange and Brutsaert 
(1989) data yield Zo/h=0.06 and do/h=0.30, while the data of Gash et al. (1986) yield 
Zo/h=0.094 and do/h=0.75. The Raupach formulation in Figure 5 passes between these 
points for Zo/h. In fact, these points fall within the scatter of the data used to construct the 
original graph [Raupach, 1994]. The Raupach formulation fordo/h on Figure 6 passes 
slightly above the Gash et al. (1986) point, but far above the Parlange and Brutsaert (1989) 
point. Also plotted on Figures 5 and 6 are the Parlange and Brutsaert (1989) data assuming 
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h=10 m. In this case, the point for zjh (i.e., Zo/h=0.12) falls somewhat higher than the 
Raupach formulation on Figure 5, but the point for do/h (i.e., do/h=0.6) is much closer to 
the Raupach formulation on Figure 6. 

A comparison of the various roughness estimates for the Landes Forest is provided in 
Table 3. The values fall generally within the broad range of published values for other 
coniferous forests (e.g. Garratt, 1977, Jarvis et al., 1975 and Thom et al.,1975) that 
indicate Zo/h to lie usually between 0.02 and 0.2, and d<>/h mostly between 0.6 and 0.9. 

Although the region analyzed consists of over three quarters pine forest, the remaining 
part contains a mixture of agriculture and shrub patches. There was no information on 
vegetation type or structure in these patches except for fractional cover, and thus, they were 
included in the analysis assuming the same overall mean canopy similarity parameter, 

N = 6.1. Results indicate that both the pine and agriculture/shrub regions possess a wide 
range of roughness values. The computed roughness of the agriculture/shrub region yields 
a mean of about 0.08 1 and a standard deviation, sza = 0.025, whereas the forested mean 
was estimated to be 0.055 and the standard deviation, szf “ 0.013. For the displacement 
height, the agriculture/shrub region possesses a mean of 0.71 and a standard deviation, sja 
= 0.121, while the forested portion possesses a mean of 0.82 and standard deviation, s<jF - 
0.061. The influence of agriculture/shrub on the roughness histograms of Figures 7 and 8 
is illustrated by drawing arrows of one standard deviation about the mean for both pine and 
agriculture/shrub. Although variance in roughness is generally not available, the computed 
means for the two vegetation types are consistent with the published values (Brutsaert, 
1982). 
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5.2 Sensitivity to fractional vegetation cover 


The sensitivity of the above formulation can be examined by rearranging Figures 5 and 
6 in terms of fractional cover. Inserting Equation (15) for the particular Landes case with 
N=6.1 directly into (3) and (5), and again applying Raupach's global values for the 
constants, yields plots of Zo/h versus m and d<yh versus m for the Landes Forest, shown in 
Figures 9 and 10, respectively. Figure 9 indicates that the greatest sensitivity of zo/h occurs 
at the lowest values of m, less than about 0. 15. At about m = 0.25, z^h reaches it peak at 
this theoretical point where the flow begins to skim over the canopy. For m greater than 
0.25, z^h is relatively insensitive to the fractional cover. 

5.3 Scale limitations 

The scales to which the frontal and canopy area parameterizations are valid are difficult 
to quantify, although the following guidelines must be met. They apply where the spatial 
distribution of the vegetation cover is statistically homogeneous and of one general 
geometric shape at the subpixel level. Thus, the upper limit of the satellite pixel scale for 
which they are applicable is variable and site specific, highly dependent on the large-scale 
spatial distribution of a given vegetation type. The lower limit depends on the relative scale 
of the pixel as compared to the scale of the individual canopy frontal area, as it is necessary 
to insure a statistically significant number of canopies within the pixel. A suitable lower 
limit criterion that meets that objective is Ap/vA t » 10. 
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6. Conclusions 


The theoretical parameterizations presented herein for the frontal and canopy area 
indices of random plants provide a mechanism for estimating local aerodynamic roughness 
over flat terrain, in terms of fractional vegetation cover observed from nadir-viewing satellite 
sensors. The formulations, graphically displayed in Figure 2 in terms of the non- 
dimensional similarity parameter, are general and theoretically can be applied to any natural 
vegetation of Poisson distribution, not only coniferous forests. The fractional cover 
estimates required in the formulations can be obtained from any source, satellite based 
approaches or other gridded or non-gridded datasets. 

When incorporated into Raupach's formulas and applied to the Landes Forest test case 
near Lubbon during HAPEX-Mobilhy, the roughness estimates compare reasonably well to 
reference values determined from field measurements taken during the experiment at the time 
of the satellite observation. They also fall within the range of roughness for coniferous 
forests. The limitations to the present analysis is that it applies to vegetation roughness, and 
the superposition of vegetation on topographic roughness is not addressed. Also, the ability 
of the parameterizations for frontal area to account for within-tree leaf density needs to be 
investigated. 

Overall, however, the good agreement between the theoretical and experimental 
estimates for the Landes Forest case warrants further testing on other natural vegetated 
systems, forested and non-forested. While satellite remote sensing can not replace direct 
field measurements, the above approach may be suitable in natural vegetated regions when 
there exists only limited ground truth. 
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Figure Captions 


Figure 1 . Schematic figure of geometric canopy model for the particular case of plants 
represented as cones on a post with Poisson distributed plant centers. Plants are 
geometrically similar. Although varying in size, they possess the same basic shape and 
constant diameter to height ratio, D/H. 

Figure 2. Canopy area index versus fractional cover for any Poisson distributed canopy. N 
is the canopy area similarity parameter, or the ratio of the total outside surface area of the 
canopy to the vertically projected area. Physically, large N represents tall plants with 
narrow canopies while small N represents short plants with wide canopies. 

Figure 3. Frequency distribution of fractional vegetation cover, Landes Forest near 
Lubbon, June 23, 1986, based on Jasinski (1996). 

Figure 4. Frequency distribution of canopy area index, Landes Forest near Lubbon, June 
23, 1986. Distribution calculated by applying Equation (15) with N=6.1 to the frequency 
distribution of fractional cover in Figure 3. 

Figure 5. Normalized aerodynamic roughness versus canopy area index. Plot compares 
Raupach's (1994) theoretical formula with estimates from field data obtained over Landes 
Forest during HAPEX-Mobilhy in June 1986. Field data are plotted versus mean canopy 
area index estimated in the present study. Data from radiosounding analysis by Parlange 
and Brutsaert (1989) are plotted using actual tree height of h - 20m as well as their 
estimated effective tree height of h = 10 m. 
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Figure 6. Normalized displacement height versus canopy area index, Landes Forest near 
Lubbon, June 1986. Plot compares Raupach's (1994) empirical formula with estimates 
from field data obtained over Landes Forest during HAPEX-Mobilhy in June 1986. Field 
data are plotted versus mean canopy area index estimated in the present study. Data from 
radiosounding analysis by Parlange and Brutsaert (1989) are plotted using actual tree height 
of h = 20m as well as their estimated effective tree height of h = 10 m. 

Figure 7. Frequency distribution of zjh, Landes Forest near Lubbon, June 1986, obtained 
by applying satellite derived canopy area index distribution of Figure 4 to Raupach's (1994) 
theoretical formulation. Arrow bars indicate one standard deviation of zjh. about its mean 
for the pine forest pixels, szf and the agricultre/shrub pixels, Sza- 

Figure 8. Frequency distribution of do/h, Landes Forest near Lubbon, June 1986, obtained 
by applying satellite derived canopy area index distribution of Figure 4 to Raupach's (1994) 
theoretical formulation. Arrow bars indicate one standard deviation of dc/h about its mean 
for the pine forest pixels, SdF and the agricultre/shrub pixels, sja- 

Figure 9. Normalized aerodymamic roughness versus fractional canopy cover. Results 
obtained in similar manner to Figure 5 except with abscissa converted to fractional cover, 
for the particular case of the Landes Forest near Lubbon, June 1986, where the mean 
canopy area index is estimated to be 6.1. 

Figure 10. Normalized displacement height versus fractional canopy cover. Results 
obtained in similar manner to Figure 6 except with abscissa converted to fractional cover, 
for the particular case of the Landes Forest near Lubbon, June 1986, where the mean 
canopy area index is estimated to be 6.1. 
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Table 1. Summary of Canopy Geometry Similarity Formulas V 
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JV Where D= maximim canopy width; H= canopy height; h s = height of post, 0 - solar zenith angle; 
X - sin-1 (D/Htan 0); p = h s (tan 0)/D; P = cos-l{[(l+(2h s /H)][(l-cos0 , )/sin0']}; 

0’ = tan-l [(Dtan0)/H]; | = cos-l[(l-cos0')/sin0']; e = [1 + (D/H)2]0.5. 

2/ Assumes D < h s tan0. 

3/ Assumes D(1 + l/cos0') < (H+2h s )tan0. 





Table 2. 

Estimated site parameters for the Landes Forest near Lubbon, June 23, 1986. 


Item (units) 


Symbol Estimate 


Solar zenith angle (degrees) 

0 

31.6 

Fractional canopy cover 

m 

0.67 

Solar-geometric similarity parameter 

T1 

1.05 

Mean canopy diameter (m) 

D 

5.1 

Mean canopy height (m) 

H 

12.8 

Mean stem height (m) 

hs 

7.2 

Mean total canopy height, H+h s (m) 

h 

20.0 

Frontal area similarity parameter 

V 

1.59 

Canopy area similarity parameter 

N 

6.1 

Mean Canopy Area Index 

A 

3.4 




Table 3. 

Comparison of roughness lengths for the Landes Forest near Lubbon, June 23, 1986. 


Source 

Height 

(m) 

zo 

(m) 

do 

(m) 

z o/h 

do/h 

Present study 

20 

1.2 

16.0 

0.061 

0.79 

Parlange and Brutsaert (1989) 

ioi/ 

1.2 

6.0 

0.12 

0.60 

II 

20 2/ 

1.2 

6.0 

0.06 

0.30 

Gashet al. (1989) 

20.3 

1.9 

15.2 

0.094 

0.75 


1/ Effective height estimated by Parlange and Brutsaert (1989). 
2/ Actual forest height. 
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